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^•-^ Abstract. In this article, we study the notion of genus expansion in the 

^ N Harish-Chandra-Itzykson-Zuber matrix model. We prove that, under suitable 

Hn hypotheses, each Taylor coefficient of the HCIZ free energy admits an TV — ^ 

^ oo asymptotic expansion in powers of N~^ whose coefficients are generating 

^«H functions for a desymmetrized version of the double Hurwitz numbers, which 

^^ we call monotone double Hurwitz numbers. We prove that the monotone 

^h double Hurwitz numbers exhibit the main structural properties of the usual 

^— ^ double Hurwitz numbers: their total generating function is a solution of the 2D 

Toda Lattice equations, and the numbers themselves are piecewise polynomial 

I I functions on pairs of partitions. 
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0. Introduction 

0.1. Overview. The Harish-Chandra-Itzykson-Zuber matrix model is a complex, 
unit-mass Borel measure /i^r on the group oi N x N unitary matrices. This measure 
is by definition absolutely continuous with respect to the Haar probability measure 
on U(A^), being given by the density 



I-n(z) 



in{z) 

where z is a complex parameter and A, B are N x N complex matrices. The 
partition function of the model, 



I^{z) = / ^^NtriAUBUn^U, 



is known as the Harish-Chandra-Itzykson-Zuber integral. This integral was first 
considered by Harish-Chandra in his study of differential operators on semisimple 
Lie algebras [5^ , where it was evaluated as a ratio of determinants under the as- 
sumption that A, B lie in the Lie algebra of the unitary group. It appears in various 
other contexts, including the study of large deviations in the spectral measure of 
Gaussian sample covariance matrices |27| , as the reproducing kernel of a distin- 
guished inner product on symmetric polynomials |49j . and especially in the theory 
of matrix models, where Harish-Chandra's formula was independently rediscovered 
by Itzykson and Zuber in the course of their study of the Hermitian two-matrix 
model [33], see also [69]. It is in this latter context that the problem of analyzing 
the A^ — >■ oo asymptotic behaviour of the free energy 

Pn{z)^ — log/w(z) 

of the HCIZ model first arose. This problem has since received considerable at- 
tention from mathematicians, notably in the work of Collins, Guionnet, and their 
collaborators, who have addressed it using the methods of large deviation theory, 
Schwinger-Dyson equations, and non-commutative differential calculus [SI HI HH HZ]- 
Apart from its importance in the asymptotic analysis of Hermitian multi-matrix 
models [231 [33] , the large N behaviour of the HCIZ free energy is of interest for 
another, deeper reason. It is well-known that the free energy of Hermitian matrix 
models admits a perturbative expansion in powers of N~'^, often called the "genus 
expansion," each order of which is a generating function enumerating polygonal 
discretizations of a compact two-dimensional surface of given topology 31 ESI US] ■ 
The genus expansion is a key ingredient in the formulation of Witten's famous 
conjecture relating two-dimensional quantum gravity to intersection theory on the 
moduli space of curves [65j. Several proofs of Witten's conjecture have been given, 
including [SnilSElHSlISS]- Theoretical physicists have postulated that the free energy 
of the HCIZ matrix model admits an analogous genus expansion with similar ties 
to algebraic combinatorics and enumerative algebraic geometry [331 ESI [55, . In this 
article, we will give a precise formulation and proof of this conjecture. 

0.2. Main results. The two prototypical ensembles of random matrix theory [U] 
are Wigner's Gaussian Unitary Ensemble (CUE) and Dyson's Circular Unitary 
Ensemble (CUE). Just as the Hermitian one-matrix model is a deformation of the 
GUE, the HCIZ model may be viewed as a deformation of the CUE. The genus 
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expansion of Bessis, Itzykson and Zuber j4] is an asymptotic expansion of the 
derivatives of the free energy of the Hermitian one-matrix model, the 5*^ coefficient 
of which counts maps on a surface of genus g. Our main result on the asymptotics 
of the HCIZ model is the analogue of the Bessis-Itzykson-Zuber genus expansion. 

Theorem 0.1. Let (An), (Bn) be two sequences of N x N normal matrices whose 
spectral radii are uniformly bounded. Suppose there exist two sequences of complex 
numbers {(f>k), (V'fc) o,nd a nonnegative integer h such that, for each fc £ N, 

as A^ — > cx). Then, for each d G N. the d*^ derivative of the free energy -Fjv(z) at 
z — admits an N -^ 00 asymptotic expansion to h terms: 



3=0 

The coefficients in this expansion are given by 



i^m-E^ + ofi^ 



a,fl\-d 



whe 



and Hg{a, (3) is the number of (r + 2) -tuples (cr, p, ti, . . . , r^) of permutations from 
the symmetric group S(d) such that 

(1) a has cycle type a, p has cycle type /?, and the ti are transpositions; 

(2) The product upri . . . r^. equals the identity permutation; 

(3) The group {a, p,Ti, . . . , Tr) acts transitively on {1, . . . , d}; 

(4) r^2g~2 + i{a)+m; 

(5) Writing Ti — (sj ij) with Si < ti, we have ti < ■ ■ ■ < t^. 



Theorem 



0.1 



is proved in Section 111 below. The convergence of Fj^'{0) was first 
established by Collins [5] , who gave a different combinatorial interpretation of the 
limit Co,d in terms of properly two-coloured planar maps, see also |69j . Theorem 0.1 
gives the full asymptotics of the Taylor coefhcients of Fn(z) under the assumption 
of rapid convergence of the moments of ^jv ^nd Bpj. For example, if there exist 
constants a,b > such that 

ltr(^^) = -0fe + O(e-'^^) 
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as N ^- cx), then Fj^'{0) admits an asymptotic expansion on the scale N~^^ to 
arbitrary order, whose combinatorial meaning is as stated above. 



The geometric meaning of the asymptotic expansion claimed in Theorem 0.1 
which justifies its designation as a "genus expansion," is best understood by pursu- 
ing the analogy with Hermitian matrix models. The genus expansion of the Hermit- 
ian one-matrix model encodes a combinatorial method for constructing a compact 
Riemann surface (or smooth projective curve), namely by glueing together polygo- 
nal tiles cut out of the complex plane. Another recipe for constructing a compact 
Riemann surface is to realize it as a branched covering of the Riemann sphere (or 
projective line) P^. Indeed, it is a consequence of the Riemann Existence Theo- 
rem that, given a finite-sheeted topological branched covering / : S* — > P^ of the 
Riemann sphere by a compact surface S, there is a unique complex structure on 
S which makes this map holomorphic. A classical construction due to Hurwitz 
[5T1 \n\ encodes a given d-sheeted branched covering as a transitive factorization of 
the identity in S(d). Such factorizations are often called constellations [40 . Hur- 
witz's construction takes as input a branched covering / : 5* — ^ P^ together with 
a labelling of the branch points of / and a labelling of the points in the fibre of / 
over a specified unbranched basepoint, and outputs a constellation whose factors 
are determined by the lifts of small loops on the sphere encircling the branch points 
of /. The cycle type of each factor coincides with the monodromy of / over the 
corresponding branch point. 

A particular case of Hurwitz's construction produces the double Hurwitz numbers 
Hg{a,f3) considered by Okounkov [17] and further studied by Goulden, Jackson, 
and Vakil [22'. Given two partitions a, /3 h d, the double Hurwitz number counts 
(up to an appropriate notion of isomorphism) degree d branched covers / : S' — ^ P^ 
of the Riemann sphere in which the source curve S has genus g and the map / 
has monodromy a over 0, /? over oo, and r = 2g — 2 + £{a) + £{(3) additional 
simple branch points at fixed positions, the number of which is determined by 
the Riemann-Hurwitz formula. Applying Hurwitz's construction, Hg{a,(3) countaj 
(r + 2)-tuples (cr, p, ti, . . . , t^) of permutations from S{d) verifying the first four of 
the combinatorial conditions listed in Theorem lO.il 

Thus, the coefficients Cg^d = '^ Hg{a, /3)(j)aipij appearing in Theorem 



0.1 



enu- 



merate certain degree d branched covers of P^ by curves of genus g. The branched 
covers being enumerated have arbitrary ramification over two given points of P^ , and 
simple branching over an appropriate number of additional fixed points. The fifth 



combinatorial condition in Theorem 0.1 which distinguishes the monotone double 
Hurwitz numbers Hg{a,f3) from Okounkov's double Hurwitz numbers, is a special 
feature of the HCIZ model whose origin will become clear below. In a sense, it is an 
elaboration of the relationship between the Catalan numbers g^-j- ( ^ ) , which count 
monotone trees, monotone parking functions, etc., and the Cayley numbers d'^~^, 
which count symmetrized structures of the same type. This analogy is developed 
in full in our second paper on this subject [17j . which contains a detailed com- 
binatorial analysis of the single monotone Hurwitz numbers Hg{a) = Hg{a, (l'')) 



proceeding in tandem with the well-developed combinatorial theory of the single 
Igia)=Hgi 



Hurwitz numbers Hg{a) = Hg{a, (I'')). 



The usual definition of the Hurwitz numbers would include a further division by d\ in order to 
compensate for reparameterizations of the domain. For our purposes, it will be more convenient 
to omit this division — this is like working with labelled maps instead of unlabelled maps. 
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Theorem |0.1| describes the N ~> oo asymptotic behaviour of the Taylor coef- 
ficients of Ff^[z) about 2 = 0. A natural question is whether this information 
describes the asymptotics of Fj^{z) itself. Indeed, this is the principal question 
in the analytic theory of matrix models, and has been intensively studied in the 
Hermitian case, see [5l[T6j[2l]- Substituting the asymptotic expansion of Fj^ (0) 
into the Maclaurin series of Fpf(z) and formally changing order of summation leads 
to the following conjectural asymptotic expansion of Fn{z). 



Conjecture 0.2. Under the hypotheses of Theorem 0.1 Fn[z) admits the N 
asymptotic expansion 



where 



' C,(z) , / 1 



o 



^23 I ]\f2h 



This asymptotic expansion holds uniformly on compact subsets of the open disc 
D(0,rcM^^), where M is the least upper bound of the spectral radii of the matrix 
sequences (An), {Bn), and re is the critical value 

2 

Conjecturc |0.2| is the HCIZ analogue of the analytical form of the genus expansion 
for the Hermitian one-matrix model established by Ercolani and McLaughlin [TB] 
and Bleher and Its [S] . Section [T] includes several steps towards a proof of this 
conjecture. First, using the results of our second paper [T7] , we prove that the genus- 



specific generating functions Cg{z) appearing in Conjecture 0.2 all have radius of 
convergence rcM~'^ , where re = 2/27 is the critical value claimed above. This is 
analogous to the Hermitian case: for example, in the Hermitian one-matrix model 
with quartic potential [H [TB] , the orders in the asymptotic expansion of the free 
energy are genus-specific generating functions for quadrangulations of a compact 
Riemann surface, and each of these generating functions has the same radius of 
convergence, namely 1/48. Second, by combining our results with those of Collins, 



Guionnet and Maurel-Segala [^, we are able to conclude that Conjecture 0.2 holds 
to first order when z is restricted to a small real neighbourhood of zero. 



Theorem 0.3. Under the hypotheses of Theorem \0.1\ and with the additional as- 
sumption that the matrices An,Bi\[ are Hermitian, there exists a non-empty open 
real interval I C [—rcM~'^,rcM~'^) such that Fn{z) converges uniformly to Co{z) 
on compact subsets of F 



are 



The techniques of [S], which yield Theorem 0.3 as a corollary of Theorem 0.1 
based on the non-commutative differential calculus developed by Guionnet and her 
collaborators in the context of Hermitian matrix models, see [351 US]- Unlike the 
partition function in the theory of Hermitian matrix models, which is a holomor- 
phic function of the deformation parameters in an open set having the origin as a 
boundary point, the HCIZ integral is an entire function of z G C. This fact suggests 
that a proof of Conjecture|G.2|based on the techniques of classical complex function 
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theory may be possible. We conclude Section [l] by sketching such an argument for 
the first order case. 

In connection with two-dimensional quantum chromodynamics, physicists have 
come to consider the enumeration of branched covers of P^ as a closed string theory 



on the Riemann sphere [23]. According to Theorem 0.1 the free energy of the 
HCIZ model approximates the partition function of such a theory under a com- 
binatorial constraint. One may study the limit object directly once it has been 
identified. Section [2] contains our second set of results, which focus on the struc- 
ture of the monotone double Hurwitz numbers Hg{a, (3). The usual double Hurwitz 
numbers enjoy several remarkable properties, the most important of which are the 
connection with integrable hierarchies of partial differential equations (171 [SI] and 
piecewise polynomial dependence on ramification type when the genus and number 
of ramification points over and oo arc held fixed (22 [Ml [SS] . The goal of Section 
[2] is to prove that the monotone double Hurwitz numbers exhibit these properties. 

Theorem 0.4. Let z, q be indeterminates, and let A — {ai, 02, . . . }, i? = {61, 62, • ■ • } 
he two alphabets of indeterminates. The formal power series 

00 d ^^^ 

H(z,g,Ai?) = E|!E'?'^ E H^{^,P)Po.{A)pp{B) 

d=l ' r=0 a,;3l-d 

is a solution of the 2D Toda lattice hierarchy in the two sets of variables pi{A) , p2{A) , . 
pi{B),p2{B),..., where 

Pi{A) = ai + 02 + ■ . ■ Pi{B) = 61 + 62 + • ■ • 
P2(A) =«i + a2 + --- P2{B)=hl + hl + ... 

are the power-sum symmetric functions in the alphabets A and B. 

This result establishes the monotone analogue of the main result of [IT]. The equa- 
tions of the Toda hierarchy yield a countable set of recurrences which uniquely 
determine the monotone double Hurwitz numbers. Furthermore, via Theorem |0.1[ 
every result about the monotone double Hurwitz numbers also furnishes informa- 
tion on the HCIZ model. Theorem |0.4| describes the evolution of the asymptotic 
behaviour of the model as the limiting moments of the matrices Ajq and Bjq vary. 
In particular, Theorem [0 . 4| generalizes a result of Zinn- Justin |67j . who proved that 
in the large N limit the free energy Fj^ becomes a solution of a limiting version of 
the 2D Toda hierarchy (the dispersionless 2D Toda hierarchy). 

Finally, we prove that the monotone double Hurwitz numbers Hg{a, f3) = 
iJg(ai, . . . , Q;m,/3i, . . • ,/3n) are piecewise polynomial functions of ai,l3j for fixed 
g, m, n. The piecewise polynomiality of the usual double Hurwitz numbers Hg{a, j3) 
was proved by Goulden, Jackson and Vakil [H] and serves to support their con- 
jecture that the double Hurwitz numbers are top intersections on some universal 
Picard variety. It is therefore desirable to formulate an analogous statement for the 
monotone double Hurwitz numbers. 

Theorem 0.5. To each triple {g,m,n) consisting of a non-negative integer g and 
positive integers m, n, there corresponds a hyperplane arrangement in M™+" and 
a collection of polynomials Pg^^^Xi, . . . , Xm, J/i, • ■ • , Vn) indexed by the chambers of 
this hyperplane arrangement such that 
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|Auta||Aut^| r} f a\ B / R a \ 
n Hg{a,P) = Pg,c(ai,---,am,/3i,---,/3n) 



|q:|! 
for all partitions a,/3 with \a\ ~ \j3\,i(a) = m,£{p) = n, and (a,/?) € c. 

We will see below that the hyperplane arrangement which determines the piecewise 
polynomiality of the monotone double Hurwitz numbers is precisely the resonance 
arrangement of j34l I57j , which determines the piecewise polynomiality of the usual 
Hurwitz numbers. 

0.3. Some remarks on our second paper. Let us close this Introduction with 
a brief discussion of our second paper on this subject |17j . 

Loosely speaking, the HCIZ matrix model is "like" the Hermitian two-matrix 
model: its free energy is a solution of the 2D Toda hierarchy that counts combina- 
torial/geometric structures with two degrees of enumerative freedom. When one of 
the matrix sequences defining the HCIZ potential, say (Bn), has degenerate limiting 
moments —ipk = ^ik, Theorem |0. 1| tells us that the HCIZ model degenerates from 
monotone double Hurwitz theory to monotone single Hurwitz theory. It should 
come as no surprise to the reader familiar with Hurwitz theory that the monotone 
single Hurwitz numbers are much more accessible than their double counterparts: 
we have explicit formulas for Hg{a) in genus g = 0, 1 piZl Theorems 0.3 and 0.4], 
and we are able to prove directly, without recourse to algebro-geometric methods, 
that an ELSV-type polynomiality property holds in all genera. 

From the point of view of random matrix theory, these results show that the 
degenerate HCIZ model is "like" the Hermitian one-matrix model: its free energy is 
a solution of the KP hierarchy, we can solve explicitly for the leading and sub- leading 
orders, and for genus two and higher we find that all Cg{z) are rational functions of 
a single explicit algebraic function of z (the analogous property for the Hermitian 
one-matrix model was conjectured in [3] and has been proved very recently by 
Ercolani [H]). Our paper [T7] gives a detailed combinatorial treatment of single 
monotone Hurwitz theory leading to the counterparts of most key results in classical 
single Hurwitz theory, thereby providing an essentially complete understanding of 
the combinatorics of the "one-sided" HCIZ model. 

In summary, these developments indicate that the combinatorics of multiplying 
transpositions subject to the monotonicity constraint underlies the theory of unitary 
matrix models in much the same way that Tutte's surgery on maps underlies the 
theory of Hermitian matrix models |26j . 
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1. The Genus Expansion 



In this Section, we present the proof of Theorem |0.1[ We begin by summa- 
rizing some basic features of the HCIZ modeh entirety of the partition function, 
A^-independent bound on the real part of the free energy, lower bound on the ra- 
dius of convergence of the Maclaurin series of Fm{z). Following this, we develop a 
unitary analogue of Wick's Lemma, i.e. a combinatorial rule for integrating poly- 
nomial functions against the Haar measure on U(iV). This derivation, which builds 
on previous work [51 [TUl SH SHI HH] j relies on the invariant theory of the unitary 
groups, and thereby leads to the combinatorics of the symmetric group, with a 
special role played by the Jucys-Murphy specialization of the algebra of symmetric 
functions. The techniques developed are used to obtain a convergent power series 
representation of the first N derivatives of Fn{z) at z = 0, from which Theorem 
O.llis obtained. 



1.1. Basic features of the HCIZ model. Since the HCIZ model is unitarily 
invariant, if we assume that the matrices A, B which define its density are normal 
then we may also assume that they are diagonal without further loss in generality. 
This assumption will be made from now on. Thus we study the sequence of complex 
functions defined by 



(1.1) 
where 



In{z) = / e'Ntr{A^UB^un^-^^ 



(1.2) 



An = 



AN) 



AN) 

''N 



Bn — 



6^^ 



AN) 



are two sequences of TV x iV complex diagonal matrices whose entries are uniformly 
bounded with least upper bound 



(1.3) 



M:=sup{|ari,...,l4"'^Ur^l,...,l4"'^|:iV>l}. 



Proposition 1.1. For any z G C, we have 

Proof. Since the Haar measure is a positive Borel measure, we have 



\In{z)\ 



^zNUiA,,UB,^U')^-lj 



< 



^zNtv(AMUBMU') 



dU < gl^NUiAj^UBi^unifi^lJ^ 



The result now follows since 



ti AnUBnU 



N N N N 

1=1 j=l 

where we have used the fact that the columns of a unitary matrix are unit vectors. 

D 



i=l 3 = 1 
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Proposition 1.2. In{z) is an entire function o/z G C, with Maclaurin series 
iNiz) = ^ NV" / {tr ANUBNU*fduY-. 

d=0 \ -^ / ■ 

Proof. We will prove that the derivative of In{z) can be computed at any point 
2: g C by differentiating under the integral sign: 

1'n{z) = [ |-e"^*^(^"^^"^*Mc/ = N /"(tr AArC/BArC/*)e"^*'(^"^^"^*M[/. 

The Taylor expansion of In{z) is obtained by repeated differentiation under the 
integral sign. 
We have 

In{z) ^ J KN{z,U)dU, 

where the kernel Kj\[{z, U) = e^^trl^JvC/SjvC/ ) jg ^^ entire function of z e C and a 
continuous function of (z, [/) e C x U(iV). Let z G C be given, and let (z„)^;^ be 
a sequence of complex numbers, each distinct from z, which converge to z. Then 
{z„} X U(iV) is a compact set and the Newton quotient 

Kn{z^,U) - Kn{z,U) 

Zn Z 

is a continuous, hence bounded, function on this set. Thus we may apply the 
dominated convergence theorem: 

^^lAzn)-I.{z)^l ^^ K.iz.,U)-K.iz,U)^^^ IIk^z^UW. 

n->-oo Zn — Z J n->oo Zn — Z J OZ 

D 



Propositions 1 1 . 1| and 1.2 together say that In[z) is an entire function of order 1. 



Consequently, by the Hadamard factorization theorem, it has the general form 



=^/4"' 



(1-4) lN{z) = e^-^X{[l-^y^'^'^ 

where cjv is a complex number depending only on N and |zj | < jzj | < ... 
are the zeros of Im{z). Let Sm C C\{z5^ , Zj , . . . } be a simply connected open 
neigbourhood of z = 0. We may define a holomorphic function Fn{z) on Sn by 

(1.5) i^^(z) = 7V-2 r :^dC, zGSn. 

Jo ^n[Q) 

The free energy F]y{z) so defined is a (scaled) branch of the logarithm of In{z) on 
Sn'- 

(1.6) e^'^~(^) = /^(z), F^(0) = 0. 
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Proposition |1.1| immediately implies an upper bound on the real part of Fn{z). 
Note that this upper bound does not depend on N. 

Proposition 1.3. For any z £ Sn, we have 

^Fn{z) <M^\z\. 

On any open disc D{0,r) C Sn, the free energy is equal to the sum of its 
Maclaurin series: 

(1.7) Fn{z)^J2Pn\0)^_, 4'^(0) = ^^ivW 

d=l ■ ^ 

The radius of convergence of this power series isrN = \z[ | , the modulus of the first 
zero of the partition function Im{z). From Proposition |1.1| and Jensen's formula 

log2,r-2A7-2 



z=0 



(see e.g. [H]), one may conclude that r^ > -^^2~^^ ^^ 

1.2. Matrix group Wick Lemma. Let G C U(A^) be a closed subgroup of the 
unitary group. A generic element of G will be denoted U — {uij). 

A function / : G — ?► C is said to be polynomial if there exists a polynomial pf in 
N'^ variables such that 

(1.8) /([/) =p/(Mii,...,MwAr) 

for all U £ G. Let A C L^(G,Haar) be the algebra of polynomial functions on G. 
This algebra admits the orthogonal decomposition 



(1.9) ^ = 0^W, 



d=0 



where A^"^' is the space of homogeneous polynomial functions of degree d. Since 



the decomposition (1.9) is orthogonal, the computation of inner products (/, 5) in 
A reduces to the case where f,g belong to the same space A^'^\ Furthermore, by 
linearity of the integral, it suffices to consider the case where /, g are monomials: 

(u»(lb(l) • • ■ Ui(^d)3id),Ui'(i)j'(i) . . . Uj'(d)j'(d)) 

(1.10) _ f 

Such inner products of monomials will be called (d + (i)-point correlation functions. 
An analogue of the usual Gaussian Wick Lemma (see e.g. [ID]) would be an 
algorithm which reduces the computation of (d + (i)-point functions to the compu- 
tation of simpler correlation functions. It is not immediately clear how to obtain 
such an algorithm, since there is apparently no analogue of the "propagator" in this 
setting. One may manufacture a propagator using the invariant theory of G. This 
gives an analogue of Wick's Lemma for G. Unfortunately, the G-Wick Lemma is 
less effective than the vector space Wick Lemma: it only reduces the computation 
of (d-l-d)-point functions to the computation of a distinguished subclass of {d + d)- 
point functions, not to (1 + l)-point functions (covariances) . Moreover, in order for 
the G-Wick Lemma to be practically useful, one must have available a complete 
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description of the G-invariants in the full mixed tensor algebra over the defining 
representation of G. Following [51 [TUl SH IMj, we will indicate how the argument 
proceeds in general, and restrict to the case G = U(A^) when necessary. 

Let V be the defining representation of G with its standard inner product, and 
let V* be the dual representation. Let {ei, . . . , cat} be an orthonormal basis of V. 
Consider the G-module Vdd = V^"^ ® (y*'^®d ^^ mixed tensors of type (d, d). 

Let P be the matrix of the orthogonal projection Vdd ~^ Vm oi the space of mixed 
tensors onto the subspacc of G-invariant tensors, with respect to the standard basis 

(1.11) ei(i) (8)---(g)ei(rf) (8)e*(i) (g) •••(g)e*(rf), i,j:[d]^[N] 

of Vdd- Since G is a compact group, the orthogonal projection Vdd — >■ Vm ™^y 
be obtained by averaging the action of G against the Haar probability measure. 
Consequently, the matrix elements of P are precisely the {d + (i)-point functions 



( |1.10[ ). Now choose a basis of the invariant subspace V^, and let A be the dim Vdd x 
dim Vj^ rectangular matrix whose columns are the coordinates of these basis vectors 



with respect to the basis ( 1.11 1. Then, since P is an orthogonal projection, it factors 
as 

(1.12) P = A{A*A)-^A* 

(this is the "outer product divided by inner product" formula for orthogonal pro- 
jections familiar from linear algebra, see e.g. |60)). 

We can extract a Wick-type formula from this factorization by taking a general 
matrix element of P on the left, and equating it with the corresponding matrix 
element on the right obtained simply from the definition of matrix multiplication. 
This will give a formula for the general {d + d)-pomt function in terms of the matrix 
elements of A and T~^, where F = A* A. Thus, in a sense, the inverse Gram matrix 
F~^ is playing the same role as that played by the propagator (or covariance matrix) 
in the Gaussian vector space Wick Lemma. 

In order for the above general reasoning to produce useful formulas, we must 
have a complete understanding of the space of tensor invariants V^. Let us now 
restrict to the case G = IJ(N), where we have an explicit basis for the space of 
tensor invariants. 

We recall that a permutation tt G S(d) is said to have a decreasing subsequence of 
length k if 7r(ii) > • • • > 7r(ifc) for some indices 1 < ii < • • • < j^ < '^^ For a survey 
of the combinatorics of increasing and decreasing subsequences in permutations, 
the reader is referred to |59j ; here we only need this definition in order to state the 
following result of Baik and Rains. 

Theorem 1.4 ([IlEO])- Let SN{d) denote the set of permutations in the symmetric 
group S(d) which have no decreasing subsequence of length A^ + 1. Then the tensors 

t:[d]^[N] 

are linearly independent and span the space of invariant tensors V^^. 
This leads to the following U(A^)-Wick Lemma. 
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Lemma 1.5 (U(iV)-Wick Lemma). For any i,j,i',j' ■ [d\ — > [N], we have 

where the sum runs over pairs of permutations a,p from the set S^id) such that 
i = i' o a, j = j' o p, and {V^^)ap is the (cr, p)-element of the inverse of the matrix 



r = 



J\l<:i<^P ^) 



(T,pSSjv(d) 

where c(7r) denotes the number of cycles in the permutation tt G S[d). 

The U(A^)-Wick Lemma was first stated in the physics literature by Samuel 
[56] : it was rediscovered and proved by Collins in [8 . The general nature of the 
argument allows it to be extended to other classical groups [10', and even to the 
setting of compact matrix quantum groups in the sense of Woronowicz, see [5]. 
These developments have led to a statistical theory of random matrices sampled 
from compact quantum groups, see 13". For our purposes, however, we will need to 
develop a finer understanding of the original V{N) case. 

Applying the U(A^)-Wick Lemma to the Maclaurin series of the HCIZ integral, 
we obtain the following. 

Proposition 1.6. The Maclaurin series of I]y{z) is 



In{z) 



d=0 ^ creSjv(d)peSjv(d) ^ 

where, for a of cycle type a, Po-(Ajv) denotes the power-sum symmetric function pa 
specialized at the eigenvalues of An (and similarly for p p{B jq) ) . 

Proof. Recall that our matrices A and B are diagonal: A = diag(ai, . . . , ap{), B = 
diag(6i, . . . , bisi), where we have omitted the dependence of these matrices on N to 
simplify notation. We have 

N N 

ti AUBU* = 2. /, ciibjUijUij, 
1=1 j=i 
and more generally 

(trAUBU*)'^ = 

E E a«(i) ■ • • (^i{d)bj{i) ■ ■ ■ bj(d)Ut{i)j(i) . . . Ui(rf)j(d)Mj(i)j(i) • • • "»(<i)j(d)- 

i:[d]^[N]j:[d]^lN] 

Applying the U(Af)-Wick Lemma to this expansion, we have 



U»(l)j(l) • ■ • U^d)3{d)U^(l)J(l) . . . Wi(d)j(d)d[/ : 



/ , / , \- JctPi 



(r-^). 



rTeSiv(d) peSjv(d) 

o-eAut(j) peAut(j) 
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where Aut(i) denotes the group of automorphisms of i under the natural action of 
S{d) on functions i : [d] ^^ [N] by permutation of arguments. Thus by Proposition 
|1.2[ we have 

In{z) 

d=0 ■ i:[d]^[N]j:[d]^[N] aeSfiid) peSfiid) 

o-sAut(i) peAut(j) 

°° z'' 
^Y^^"" Y Y (r~')ap Y Y ^^{l)■■■a^id)bJ(,l)■■■bJ(^d) 

d=0 ■ CTeSjv(d)peS„(d) i:[d]^[N]j:[d]^lN] 

j:eFix((T) jGFix(p) 

- ^d 



E;7[^' E Y ir-\pPAA)Pp{B), 



d\ 

wherepo.(A),pp(i3) denote the power-sum symmetric functions labelled by the cycle 
types of a and p specialized at the eigenvalues of A and B. 

D 

1.3. Resistance matrices. The propagator in the U(iV)-Wick Lemma, i.e. the 
Gram matrix associated to the orthogonal projection of Vdd onto its subspace of 
U(A'^)-invariants, may be interpreted as a specialization of the resistance matrix 
associated to a certain metric on the symmetric group S{d). 

Let T be a set of transpositions which generates S(d), and let |7r| = \'k\t be the 
corresponding word norm on S{d). Thus \cFp~^ \ is the length of a geodesic joining p 
to a in the Cayley graph of S(d) corresponding to the generating set T. This defines 
a metric on S(ci). 

Let (7 be a complex variable. The resistance matrix 51 = ^riq) associated to the 
T- induced metric on S{d) is by definition the d\ x d\ matrix 



(1.13) n 



qWp- 



a,peS{d) 



If we think of the Cayley graph as an electrical network in which each edge is a 
wire of resistivity q, then the (a, p)-entry of f2 represents the compound resistance 
encountered by a charge traversing any geodesic from p to a. 

One may inquire after various properties of the resistance matrix associated to 
a given metric on S(d), such as its determinant as a function of q, or the form of 
the inverse matrix fl~^ away from the zeros of the determinant. Questions of this 
sort were posed by Zagier [55] in the case that T — {(1 2), (2 3), . . . , (d — 1 d)} are 
the Coxeter generators of the symmetric group. 

For the Coxeter metric, \tt\ is equal to the number of inversions in tt. Zagier used 
this fact to reduce the problem of proving the existence of a Fock space represen- 
tation for the g-Heisenberg algebra when — 1 < g < 1 to the problem of proving 
that the resistance matrix fi is non-singular in this range. He then proved the 
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invertibility of Q away from certain roots of unity by observing that it is the image 
of the resistance element 

(1.14) ^= Yl «'''''' 

Ti-GS(d) 

in the regular representation of the group algebra C[S(rf)], with respect to the 
standard (permutation) basis. This provides leverage on the problem since, as 
shown by Zagier, the resistance element can be factored in C[S((i)] as 

(1.15) t^ = CiC2-.-Cd, 

where the factors Q are certain group algebra elements whose images in the irre- 
ducible representations V^ of C[S((i)] are easy to understand. Knowledge of the 
action of the individual factors Q in irreducible representations can then be as- 
sembled into knowledge of the action of a; in the regular representation using the 
isotypic decomposition 



(1.16) C[S(d)] = 0(dimA)l^^ 



Ahd 



of the group algebra. Indeed, this is precisely the strategy invented by Frobenius 
in his solution of Dedekind's group determinant question, which was the initial 
motivation behind the character theory of finite groups, see [35] as well as the 
discussion in [66] . In this way Zagier obtained the beautiful formula 



d-l 



(1.17) detn^Hil-q^^^+'^r, e.= M J(*-l)!(d-z)!, 



4 = 1 



from which it is clear that fl is non-singular provided q is not an i{i + l)-st root of 
unity for any i = 1, . . . ,d ~ I. 

Later, it was pointed out by Hanlon and Stanley |28j that the resistance matrix 
fl may be viewed as a specialization of the Varchenko matrix associated to the Ad-i 



hyperplane arrangement, and so the determinant evaluation (1.171 also follows from 
general formulas due to Varchenko [64] . 

In our setting, we are interested in the resistance matrix corresponding to the 
generating set T = C(2id-2-), the entire conjugacy class of transpositions. The 
corresponding word norm is given by |7r| = d — c{tt), with c(7r) the number of cycles 
in TT (including 1-cycles). Thus, for d < N, the propagator in the U(A^)-Wick 
Lemma is given by 

(1.18) r = N'^n, 

with resistivity specialized at q = 1/N. 

In order to understand the propagator in the U(A^)-Wick formula, we may follow 
Zagier's strategy and factor the corresponding resistance element w € C[S((i)] into 
tractable pieces. The factorization implementing this strategy has been known since 
the work of Jucys [3S], see [131 Proposition 2.1] for a proof. 
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Lemma 1.7. For the all-transpositions distance on S(ci), the resistance element ui 
factors as 

a; = (1 + q.h){l + gJs) . . . (1 + qJd), 
where Ji :— and 

Jt = j2{st) 

s<t 

for2<t<d. 

The transposition sums appearing in the above lemma are the J ucys- Murphy 
elements: 

J2 = (1 2) 
J3 = (1 3) + (2 3) 
(1-19) J4 = (14) + (24) + (34) 



They can ahernatively be written as 

(1.20) Jt = y^ transpositions in S(i) — \^ transpositions in S{t — 1), 

which shows that Jt belongs to the (maximal) commutative subalgebra of C[S((i)] 
simultaneously generated by the images of the class algebras Z{1), -2(2), . . . , Z{d) 
under the standard embedding S(fc) ^^ S((i) for A: < d (this is the Gelfand-Zetlin 
algbera; see [55]). Note also that the right hand side of the factorization in Lemma 
|1.7|is precisely the generating function for the elementary symmetric functions, 



(1.21) Yl l''^k{xi,X2, ...) = 1[{1 + qxi), 

r=0 4=1 

specialized on the commutative alphabet S^ = {{ Ji, J2, . . . , Jd, 0, 0, . . . }}. Thus an 
equivalent statement of Lemma |1.7| is 



(1.22) efe(Sd)= Y. C-M, 

^\-d 
e(p)=d-k 

where C^ € Z{d) is the sum of all permutations of cycle type fi. It follows from 
the fundamental theorem in symmetric function theory, A = C[ei, 62, . . .], that 
/ n- /(Sd) defines a specialization A — > Z{d) mapping the algebra of symmetric 
functions onto the class algebra. 

Since /(S^) S Z{d) for any symmetric function /, it follows from Schur's Lemma 
that /(5d) acts as a scalar operator in any irreducible representation V^ oiC[S{d)]. 
It is a remarkable result of Jucys [3S] that the central character of /(S^) acting in 
F'*' is given by the substitution rule /(S^) i~> f{A\), where 

(1.23) A^ = {{c(n) : D e A}} 
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is the alphabet of contenttjjof the Young diagram X\- d. In particular, the central 
character of the resistance element uj in V"^ is 



(1.24) Hil + qciD)). 

We now see from the isotypic decomposition C[S((i)] — 0^|_^(dim A)!/^"^ that the 
determinant of the resistance matrix Q corresponding to the all-transpositions word 
metric on S(d) is 



d-l 

(1.25) detn=Y[{l-c'^q^y'', ec=^dimA. 

c=l Xhd 

ceAx 

It follows from the above computation that the resistance matrix fl associated 
to the all-transpositions distance on S{d) is invertible for q ^ {±l,±i, . . . ,±-£-^}. 
The entries of the inverse matrix fl~^ are rational functions of q which are analytic 
in the disc \q\ < ^ztj;. In particular, by virtue of the reciprocal relationship between 
the elementary symmetric functions e^ and the complete symmetric functions 



(1.26) hr{xi,X2,-.-) = ^ Xtj^Xt2---Xt^, 

ti<t2<---<tr 

the Maclaurin series of the matrix elements oi fl^^ is 



(1.27) (f7-i).p = ^(-9)^qap-i]/i,(S,) = (-l)l-''"'l^g'-[ap-i]/i,(Srf), 

r=0 r=0 

with radius of convergence at least g— j- for each matrix element. Note that the 
above is not an alternating series; its coefficients are either non-negative integers 
or non-positive integers, depending on the parity of the permutation ap~^. 

If we interpret fl~^ as the "propagator" in the U(iV)-Wick formula, we see a ma- 
jor departure from the Gaussian Wick Lemma: the propagator makes power series 
contributions. This feature of polynomial integrals over V{N) was first pointed out 
by De Wit and 't Hooft [12 . The poles of the matrix elements of il^^ are known 
as "De Wit - 't Hooft anomalies" in the physics literature [55] , 



1.4. Asymptotic expansion of Taylor coefficients. We now substitute the 
power series form of the propagator il^^ = N'^T^^ into the series expansion of 



the HCIZ integral obtained in Proposition 1.6 This yields: 



Recall that the content of the cell D in row i and column j of a Young diagram Aisc(n) = j- 
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In{z) = X! ( X! X! (^~^)'^P-P^(^)Pp(^) ) "jj" + higher terms in z 
= 11 \Y^ HH\-^] {^P~^\K{^d)vM)Vp{B) \ — + higher terms in z 

d=0 ^ aeS{d) peS{d) r=0 ^ ^ ' 

N , oo / 1 \ "^ \ <i 

= H,\ H, X! ( ~ 77 ) [C(id)]C„C^/ir(Srf)pc.(^)p/3(-B) ] -^ + higher terms in z. 



We have arrived at the identity 



In{z) = 12[ 12 H \n) i'^{i'')]'^''^^l3^r{'^d)Pa{A)pp{B) j — +higher terms in z, 

which gives an absolutely convergent series expansion for each of the first N coef- 
ficients in the Maclaurin series of In{z) in terms of multiplication in the center of 
the symmetric group algebra. 

From the definition of the JM elements and the complete symmetric function hr, 
one sees that the expression [C(^id-j]CaCphr{'B,d) counts (r+2)-tuples (cr, p, ri, . . . , r^) 
of permutations from the symmetric group S(d) such that 

(1) (7 has cycle type a, p has cycle type /3, and ti, . . . , r^ are transpositions; 

(2) The product apri . . .t^ is the identity permutation; 

(3) The transpositions ri — {si ti), . . . ^t^ — {sr tr) satisfy ii < • • • < tr. 

Observe that z is an exponential marker for the size d of the ground set in this 
combinatorial problem, while —l/N is an ordinary marker for the number r of 
transposition factors, which must be ordered according to the monotonicity con- 
straint ti < ■ ■ ■ < tj.. Since 

(1.28) e^'^~(-) = /Ar(z), Fjv(0)-0, 

it follows from the general theory of generating functions [THl Chapter 3] that the 
Maclaurin series of Fpf{z) is 

Fn{z)=N-^Y.[ E E(-]^) ^'■(a,/?K(^)pMi?))|j+higher terms in z. 

where H^{a,P) is the number of solutions to the same combinatorial problem as 
above, but now with the additional condition that the subgroup of S(d) generated by 
the factors cr, p, ri, . . . , r^ acts transitively on the points {1, . . . , d}. This is exactly 
like the passage from all maps to connected maps via the logarithm in Hermitian 
matrix models. By the Riemann-Hurwitz formula, -ff'"(a,/3) = Hg{a,l3) when 
r = 2g — 2 + i{a) + £{f3) for g > 0, and vanishes otherwise. We thus obtain the 
following convergent series representation for each of the first N derivatives oi F^{z) 
at z = 0. 
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Theorem 1.8. For I < d < N, we have the absolutely convergent series represen- 
tation 



where 



N^9 ' 



a,l3\-d 



Theorem [OTT] may now be deduced from Thcorem |1.8| as follows. Our assumptions 
on the uniform boundedness and convergence to limiting moments at a specified 
rate of the matrix sequences (A^), (Bn) are equivalent, respectively, to the uniform 
bounds 

(1 29) 

and the N ^ oo asymptotics 

(1.30) ) 

for all partitions a, P and a specified /i G {0, 1, 2, . . . } U {oo}, where 

V'/3 = n ^ft 



and taking h = oo should be interpreted as saying that the equations ( 1.30 ) hold for 
any non-negative integer h, which would be the case, for example, if the moments of 
A^ and Bn tend to their limits exponentially fast. In particular, the asymptotics 



(1.30) say that the quantities Cg^d.N arising in the absolutely con verge nt series 
representation of the Taylor coefficients F^ (0) obtained in Theorem 1.8 satisfy 



(1-32) (^g.d.N = Cg^d + o 



I JY2/i 



Prom 1.8 we have 



N"^ (f^'^ (0) - X: %f ) = C,.d^N - ^'^±i^ 



3=0 ^ 



for all N > d. Now 
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a,l3\-d 

a, ,3 
< M2'^p(d)2(d!)29-2+2'i, 

where p{d) denotes the number of partitions of d and we have used the estimate 

which fohows immediately from the definition of the monotone double Hurwitz 
numbers. Furthermore, we have 

Ck+i4,n . CU+24.N . ^ APMd)'id\r^'+'^ ( jdiy {d\y \ 

N^ m + • • • - Ar2 \ N^ Ni J 

M2''p(d)2(d!)2(''+'') 1 



N2 1 _ Mi 



for all A^ > d!, so that 



^ 9=0 ^ \ 

as A^ — > oo. The result now follows from the fact that 

(1-33) Cg^d.N ^Cg^d + oij^ 

as A^ — > oo. 

1.5. Convergence of Hurwitz generating functions. Recall that the free en- 
ergy of the HCIZ model is equal to the sum of its Maclaurin series in a neighbour- 
hood of z = 0: we have 



(1.34) F^{z) = yFlf\0f 



d 



d=l 



r'^ ' \Tir\rtaiVCi -y^ •' 



for \z\ < rN — \zi I, where zl ' is the first zero of the HCIZ integral. Replacing 
the Taylor coefficient Fj^ (0) with its asymptotic expansion and working formally 
with the resulting expression, 



°° ,,. yd °° / °° r .\ yd. oo . oo j 

(1.35) M.)=E^?(»)s-i: E^ s-i:]^ i:c'4 

d=l d=l ^ g=0 ' g=0 ^ d=l 



leads to the formulation of conjecture |0.2[ As a first step towards validating this 
conjecture, one must prove that the power series 
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(1-36) Cgiz)^J2^9/^' 9>0, 

d=l 

all converge on the open disc D{0, TcM^"^). We will prove this here. 
By the boundedness of spectral radii hypothesis, we have the estimate 



^ Hg{a,fi)ct^^^p <M^^ Y, Hgia,(3), 



a.fShd 



(1-37) |Cg,rf| = 

a,l3\-d 

so that our convergence problem reduces to locating the dominant singularities of 
the genus-specific generating functions 

oo ^ \ d 

(1.38) n,iz)^J2( E ^s("'^)^ 

of the monotone double Hurwitz numbers. The singularity analysis of fixed-genus 
generating functions of Hurwitz numbers is an interesting problem which seems not 
to have been addressed in the literature on Hurwitz theory. We will solve it by 
reducing to the case of simple Hurwitz numbers, and then applying a specialization 
of the results in our second paper ^/Tj . 

Definition 1.9. The monotone simple Hurwitz number Hg.d := Hg{{l'^),{l'^)) is 
equal to the number of r-tuples (ri, . . . , r^) of transpositions from the symmetric 
group S(d) such that: 

(1) The product ti . . . r^ equals the identity permutation; 

(2) The group (ri, . . . , r^) acts transitively on {1, . . . , d}; 

(3) r^2g-2 + 2d; 

(4) Writing Tj — (s^ tj) with Si < ti, we have ti < ■ ■ ■ < tr- 

The monotone simple Hurwitz number Hgd counts, up to isomorphism, a combi- 
natorially restricted subclass of the set of genus g, degree d branched covers of P^ 
with simple ramification over r = 2g — 2 + 2d fixed points of the sphere and no 
other branching, where the number r is determined by the Riemann-Hurwitz for- 
mula. The fixed-genus generating function encoding the monotone simple Hurwitz 
numbers in genus g is 

°° z''- 

(1-39) H^(z)=5]i?,,d-. 

Theorem 1.10. The generating functions Hg(z) andii^Az) have the same radius 
of convergence. 

Proof. Clearly, we have that 

^ Hg{a,f3)>Hg^d, 

a,l3\-d 

so the radius of convergence of ilg{z) is at most the radius of convergence of Hg(z). 
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Conversely, a straightforward combinatorial argument involving successive mul- 
tiplication by appropriately ordered cut transpositions I17j leads to the inequality 

iJg(a,/3) <Hg{2l'^-^,2l'^-^) for alia, /3hd 



as well as the identity 



so that we have 



Hg{21^-^,2l''-^) = ^Hg^d, 



a,p\-d 



where p{d) is the number of partitions of d. By the Hardy- Ramanuj an formula, 

Pid) = ' r- +o(l), d^oo, 
we have 

d-foo \^ 4 

Thus the radius of convergence of Hg(z) is at least the radius of convergence of 

D 



We can now combine Theorem 1.10 with the results of our second paper [T7] to 
determine the dominant singularities of the generating functions Hg(z). 

In genus (7 = 0, this is a very direct and tangible calculation. Specializing 
a — (f*) in [17, Theorem 0.3] yields the exact formula 

(1.40) ^"'"^ ^'"' ^^^-^ 



d\ d^{2d-l)\d-l 

for the genus zero monotone simple Hurwitz numbers, from which we directly find 
that the radius of convergence of Hg(z), and hence Ho(2;), is 2/27. Note that this 
value was also determined by Zinn- Justin |67) using the dispersionless Toda formal- 
ism. 

In genus 5 > 1, specializing ^7\ Theorem 0.5] at a = (1'') yields the following 
rational form. 



Theorem 1.11. Let 

I ou z, , 

77, V 77, — 1 / 

n=l 

be the unique formal power series solution of the functional equation 



„— 1 \ / 



s = z{l-2sy^ 
in the ring C[[z]], obtained by Lagrange Inversion. 
For g — 1, we have 
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H?(z) = ^log(l-2s)-^log(l-65), 
while for any g > 2, we have 

1 ^ c 6'^("^s^("^ 

H^(^) - -cs,(0) + (1 _ g,)2 E E (T- 6s)^(") ^ 

^ ' d=Oa\-d. ^ ' 

where Cga € Q are rational constants. 



It is clear from the genus zero formula (1.401 and Theorem 1 1 . 1 1| that each gener- 



ating function H^(z), g > 1, has radius of convergence equal to that of s{z), namely 
2/27. It thus follows from Pringsheim's theorerrrl that the point r^ — 2/27 is the 
common dominant singularity of all the generating functions Hg(z), 5 > 0. 

Remark 1.12. Considering in more detail the genus zero case, Zinn- Justin [67[ p. 
425] observes that by comparing with results in Tutte [63] , one finds 

(1,41) 4hs(.) . . + ,p(,) ^ I ^^ (^;:f ).^ 

where p{z) is the generating series for rooted planar maps with respect to vertices. 
Tutte specifies p{z) parametrically by 

(1.42) p = 0(l-20), (t>^z{l + 24>f. 

The results in our second paper il?, Section 3] specify Hg(z) parametrically by 

(1.43) Uz^^ - l\z^ni{z) ^s^s\ s = z(l - 2s)-2. 
These parameterizations can be reconciled by the transformation 

1.6. Convergence of the free energy. We now know that the generating func- 
tions 

(1-45) C,{z) = Y.^,/-, 

d=l 

which encode the orders of the N ^ 00 asymptotic expansion of the Taylor coeffi- 
cients of Fm{z) about 2 = 0, all converge absolutely in the open disc -D(0, TcAf^^), 



where re is the critical value 2/27. Theorem 0.3 now follows immediately by combin- 
ing Theorem |0.1| with the work of Collins, Guionnet and Maurel-Segala [9, Theorem 
0.1], which establishes the convergence of a very general class of unitary matrix in- 
tegrals in a small real neighbourhood of z = 0. We now present an argument which 
proves the convergence of the free energy Fn(z) to the generating function Cq{z) on 



Pringsheim's theorem: A power series J^ c^z"* which has non-negative real coefficients and 
radius of convergence r necessarily has a singularity aX z = r. 
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a small complex neighbourhood of the origin, modulo a technical condition which 
we have been unable to rigorously verify. 



Recall that the free energy F^{z) as defined by (1.5) is holomorphic on the 
open disc -D(0, rj\i), where r^ = \z[ | is the modulus of the first zero of the HCIZ 
integral. 

Conjecture 1.13. Zero is not a limit point of the sequence (r^r). 

Let Z? C C be an open set in the complex plane, and let Hol(D) denote the 
algebra of holomorphic functions f : D ^ C We recall that a set J-" C Hol(I?) is 
said to be a normal family if, for any compact set K G D, there exists a constant 
Mk > such that 

(1.46) II/IIk < Mk for all / e J", 

where 



(1-47) ||/|U = sup|/(z)| 

zeK 

is the sup-norm of / on K. 

Proposition 1.14. // Conjecture \1-13\ holds, then there exists < r < re and 
A'o e N such that {Fn ■ N > Nq} C Ho1D(0, r) is a normal family. 

Proof. Conjecture |1.13| is equivalent to the existence of a real number < r < re 
and an integer Nq ^ N such that Fn is holomorphic on the open disc £>(0, r) for all 
N > Nq. It remains to verify that the family {Ff^ : N > Nq} is uniformly bounded 



on compact subsets of D(0,r). This in turn can be deduced from Proposition 1.3 
and an application of the Borel-Caratheodory inequalitjQ as follows. 

Let < ri < r2 < r, and let z £ D{0,ri) be arbitrary. Then, for N > Nq, we 
have 

\Fn{z)\ < ^*^| , sup ^Fn{z) < ^^^ sup 3?Fjv(z) < "^Mll^^ 
r2 - \z\ 1^1=^2 f2 - ri 1^1=^2 r2 - ri 



where the first inequality is Borel-Caratheodory and the last is Proposition 1.3 
Consequently, the family {Fm '■ N > Nq} is uniformly bounded on the closed disc 
D{0,ri), as required. D 

Theorem 1.15. // Conjecture \l . 1 3\ holds, then fjv converges uniformly to Cq on 
compact subsets of D{0,r). 



Proof. By Proposition 1.14 and Vitali's Theorenj^ it suffices to prove that F[<, 



converges pointwise to Cq on D(0,r). 

Let zq £ D{0,r) and e > be given. We must prove that there exists N^ £ N 
such that 

N>N, =^ |Fjv(zo)~Co(zo)l<e. 



This classical inequality bounds the modulus of a holomorphic function on the interior of a 
disc by the supremum of its real part on the boundary circle, see e.g. |41j . 

Vitali's theorem: A normal family {-Fjv} C }io\(D) which converges pointwise on D con- 
verges uniformly on compact subsets of D. 
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Let |zo| < ri < r. Assuming that Conjecture 1.13 holds, we have 



E 



Id 



d=l ' d=E+l 

for any N > No and all E E N. Let < ri < r. Then, by Cauchy's estimates, we 
have 

'^ \r('^hn\ n 1^ \\FN~Co\\r, ^ ||Fjv||,.^ + ||Co||ri 



d! 



|FJV"^(0)-Co,rf|< "^" , -"'^ < 



r 



d 



for all iV > A^o and d e N, where 



llFjvlln = sup |FAr(2:)| 
|z|=ri 

is the sup- norm of Fj^ on the circle \z\ = ri, and similarly for ||Co||ri • In particular, 
our estimate on \Fisi{zq) — Cq{zq)\ becomes 



d\ 1 _ M 



N>No ^=^ |FAr(zo)-Co(zo)| < Y. 1^^ (O)-^o 

d=l "' "^ ^ 

valid for all N > Nq and any ii^ € N. Let ri < r2 < r. Applying the Borel 
Caratheodory inquality as in the proof of Proposition |1 .141 we have 

N>No => WFnWt, < -^^^^— sup ^Fn{z) < '^Mn^^ 
Returning to our estimate of [-^^(^o) — Co{zo)\, we thus have 



A^ > Ao =^ \FNizo)~Coizo)\ < J2 l^w (0)-C, 



0,d| 



Mo\\T^ + \\Co\\r.f\zo\ 



dl 1 _ M \ ri 

d=l -^ ri ^ ■ 

valid for all E E N. Now, there exists Eq G N depending on |zo|, Af, e, ri, ||Co||ri 
and r2, such that 



E+l 



2Mrir2 , ||/^_| /i I \ -E+l 



f>eo=. -- ^,!!.f°"-- r^) <./2. 



Thus we have 



2^ _ kol \ n 



A^>Ao ^ |i^Ar(zo)-Co(zo)|<X;i4'^(0)-Co,d|^+£/2. 



d\ 

d=l 



Now limAf^oo |i^lf^ (0) - Co,rf| = for all d e N, by Theorem 
A^i depending on |zo|j£j and Eq such that 



0.1 



Hence there exists 



Eq I \d 

N>m ^Y. \Pn\Q) ~ ComV-^ < e/2. 

d=l 

Setting N^ = maxjA^o, A^i} completes the proof of the pointwise convergence of Fn 
to Co on D{0,r). D 
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2. Structure of Monotone Double Hurwitz Numbers 

Having established that the monotone double Hurwitz numbers arc the combi- 
natorial/geometric objects underlying the notion of genus expansion in the HCIZ 
model, it is of interest to gather as much information as possible regarding the 
structure of these objects. In particular, it is natural to seek analogues of the 
structural properties of the usual double Hurwitz numbers in this new setting. We 
are thus motivated to search for integrable properties of the Witten-type generating 
function for monotone double Hurwitz numbers in all degrees and genera, and to 
look for polynomial behaviour exhibited by the monotone double Hurwitz numbers 
themselves, regarded as functions on pairs of partitions. 

2.1. Integrable hierarchies. 

2.1.1. The link between Hurwitz theory and integrable systems was suggested by 
Pandharipande [5T, who showed that the (at the time conjectural) Toda equation 
for the Gromov-Witten potential of P^ implies a Toda equation satisfied by the 
generating function 

(2-1) H«(-''?) = E^E77^^^ 

of the classical simple Hurwitz numbers in all degrees and genera. Here iJ^' = Hgd 
counts branched covers of P^ by curves of genus g with r = 2g — 2 + 2d simple 
ramification points at fixed positions. Pandharipande's conjecture was settled in 
short order by Okounkov 47 , who proved the considerably stronger result that the 
richer generating function 

OO ^ OO J. 

(2.2) H(z, 9, ^, i?) = E ^ E 7! E ^''("' P)Po.{A)pp{B) 

d=l ' r=0 ' a,p 

for the double Hurwitz numbers H^{a, (3) — Hg{a, I3),r — 2g — 2 + £{a) + i{(3), in 
which arbitrary ramification is permitted over and oo in addition to the r simple 
ramification points, satisfies the entire 2D Toda lattice hierarchy of Takasaki and 
Ueno lET] in the variables pi{A),p2{A), . . . ,pi{B),p2{B), . . . , where 

Pi{A)=ai + a2 + ... pi{B)=bi+b2 + ... 
P2{A)=al + al + ... p2{B)=bl + b?. + ... 



are the power-sum symmetric functions in auxiliary sets of variables A — {oi, a2, . . . } 
and B = {6i, &2, ■ • • }■ This evidence of structure in the double Hurwitz numbers 
is what led to their further investigation by Goulden, Jackson and Vakil [22] and 
many other authors since, see the discussion of piecewise polynomiality below. 

In this section we prove Theorem |0.4[ which is the monotone analogue of Ok- 
ounkov's result cited above. Introduce the generating function 

OO rf OO 

(2.3) f{z^q.A,B)^Y.^Y.l'' E [^(i'')]^"^^'^-(-'i)P"(^)P^(^)' 

d=0 ' r=0 a,0\-d 
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which is the partition function counting all (i.e. including possibly disconnected) 
monotone double Hurwitz covers of P^ . As in Section [T] the coefficient 



(2.4) 



-q^p^iA)pp{B) 



'{z,q,A,B) = [C(id)]C„C/3/i^(Sd) 



of this formal power series is equal to the number of (r + 2)-tuple3j (cr, p, ri , . . . , r^) 
of permutations from the symmetric group S(d) such that: 

(1) (7 has cycle type a, p has cycle type /3, and the Ti are transpositions; 

(2) The product apri . . .Tr equals the identity permutation; 

(3) Writing r^ = (s^ t^) with Si < ti, we have ti < ■ ■ ■ < t^- 

The variable z is an exponential marker for the size d of the ground set {1, . . . , d}, 
while q is an ordinary marker for the number of transposition factors, which must 
be ordered as in condition (3) above. Note that this differs from the generating 
function for possibly disconnected classical Hurwitz numbers, which is exponential 
in both the degree of the covering and the number of simple ramification points. 
By the Exponential Formula (a.k.a. Moment-Cumulant Formula), we have 

(2.5) \ogT{z,q,A,B)=U{z,q,A,B)=Y,^_Y.'i'^ E H^{a,(3)p^{A)pp{B). 

d=l ■ r=0 a,l3\-d 

As in [17], we will give a representation-theoretic proof that f{z, q, A, B) is a tau 
function of the 2D Toda lattice hierarchy in the variables pi{A),p2{A), . . . ,pi{B), 
P2{B), .... The basic ingredient required in the argument, namely Jucys' result on 
the spectra of symmetric functions of JM elements in irreducible representations 
of the symmetric group, was already introduced and utilized in Section [I] where 
we proved that the asymptotic expansion of the d-th Taylor coefficient Fj^'{0) of 
the HCIZ free energy is essentially the ordinary generating function of monotone 
double Hurwitz numbers in fixed degree d. 

Let 

(2.6) ^^(^^^l^'^'d^ 

denote the central character of the conjugacy class C^ G -2(d) acting in V^, and 
recall that the eigenvalue of hr{Sd) acting in V^ is hr{A\). From the isotypic 
decomposition of C[S((i)] we have 



(2.7) 



Xhd ^ 



d\ d\ 

\\-d 



Using this fact together with the change of basis from the Schur functions to the 
Newton power-sums. 



We hope that the use of the Greek letter t for both tau function and transposition will cause 
no confusion. 
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(2.8) sa(a) = ^^x;:pp(a), 

fi\-d 

we see that the generating function r may be rewritten as 



(2.9) r(z, q,AB)^J2[I[ i ac(n) •^a(^)^^(^)' 

where the summation is over ah partitions A and we have used the generating 
function 



(2.10) y2Kix,,x2,...)q'-^r\y^ 

for the complete homogeneous symmetric functions. 



The formula ( 2.9 1 already suffices to identify the formal power series t(z, g, A, i?) 
as a tau function of the KP hierarchy. In particular, it is known that a formal power 
series of the form 

(2.11) r = ^n.SA(^) 

A 

is a tau function of the KP hierarchy in the variables pi(j4),p2(^), ... if and only 
if the coefficients 1a satisfy the Pliicker relations, see e.g. [11]. The coefficients Y\ 
are called the Pliicker coordinates of the corresponding tau function. The Schur 
functions themselves satisfy the Pliicker relations and, as shown in [20] so do the 
products 



(2.12) Yy^={Wy,(U)]sx[B), 

^DeA 



where {yoi2/±ij2/±2, . . . } is an auxiliary set of variables. According to (2.9 1, the 
generating function r(z, g, A, B) is precisely of this form with 

(2.13) yc(n) 



1 - qc{U) ■ 



Thus H(z, g, A, B) = log t(z, g, A^ B) is a solution of the KP hierarchy in the vari- 
ables pi( A), p2(^), ... . 

We now wish to prove the stronger statement claimed in Theorem |0.4[ namely 
that H(z, g, A^ B) is a solution of the 2D Toda lattice equations. Let us recall that 
a formal power series solution of the 2D Toda hierarchy in the two sets of variables 
Pi{A),p2{A), . . . ,pi(S),p2(S), ... of the form 



(2.14) log(^yASA(A)sA(B 
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is referred to as a diagonal solution. The following recent result due to Carrell [6] 
asserts that diagonal series whose coefficients Yx are shifted content products are 
solutions of the 2D Toda hierarchy. 

Theorem 2.1 ([6). The sequence of formal power series 



Tn{z,q,A,B) = > Yx{n)sx{A)sxiB), n e Z, 



X 



whe 



r„(A) = e„ Y[ Vn+cio), n G z, 

1, 71 = 

is a sequence of tau functions of the 2D Toda lattice equations in the variables 
Pi{A),p2{A),...,p,iB),p2iB),.... 



In light of Theorem 2.1 the same argument we used to demonstrate that H(z, q, A, B) 
is a solution of the KP hierarchy in the variables pi {A) , p2 (A) , . . . also implies that 
H(z, q, A, B) is the n = term of a sequence of diagonal solutions of the 2D Toda 
lattice equations in pi{A),p2{A), . . . ,pi{B),p2{B), .... 

2.2. Piecewise polynomiality. The remarkable ELSV formula [TJ, which ex- 
presses the single Hurwitz numbers Hg{a) as integrals over the Deligne-Mumford 
compactification M.g^rn of the moduli space of genus g curves with m marked points, 
implies the existence of a family of polynomials Pg{xi, . . . , Xm) such that 

(2.15) \^^Hg{a)^C{g,a)Pg{a^,...,a^) 

\a\\ 

for all partitions a with £{a) — m, where the combinatorial prefactor C{g, a) is 
given explicitly by 

(2.16) C{g,a)^{2g-2 + m + d)\l[^. 

Goulden, Jackson and Vakil |2'l"i proved that this polynomiality property of single 
Hurwitz numbers is equivalent, for genus g > 2, to the rationality of the generating 
function 

.2^/^ H,(a) ^,,,\z' 



(.17) H.,,., = ^^^^^^^^_^^^^^^^,,(.,J- 

in terms of an implicitly defined set of variables obtained by Lagrange inversion. 
This rational form had been conjectured by Goulden and Jackson |19j (and proved 
in genus (7 = 2,3) prior to the advent of the ELSV formula, but continues to resist 
a proof independent of ELSV. 
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In the sequel to this paper, we prove by direct combinatorial methods ^17i The- 
orem 0.5] that the generating function 

(2.18) H,(z, A) = ^ lY^H,{a)p^iA)]- 

for the single monotone Hurwitz numbers in genus 5 > 2 is rational in an implicitly 
defined set of Lagrangian variables. This is equivalent to an ELSV-type polynomi- 
ality property for monotone single Hurwitz numbers ^17, Theorem 0.6], obtained 
without recourse to an analogue of the ELSV formula: there exist polynomials 

Pg(a;i, . . . , Xm) such that 

(2.19) ^^ij^(«)^c(a)P^(ai,...,a™), 

for all partitions a with Hia) — m, where the genus- independent combinatorial 
prefactor is given explicitly by a product of central binomial coefficients 

(2.20, C(o, = ft (t) 

over the parts of a. 

Polynomiality does not persist for the double Hurwitz numbers Hg{a,P), whose 
structure is much more complicated than that of the single Hurwitz numbers. How- 
ever, Goulden, Jackson and Vakil [52] showed there is a suitable replacement for 
polynomiality in this context: piecewise polynomiality. The flavour of this result is 
as follows. For fixed m,n we may view pairs of partitions (a,/3) with \a\ = |/3| and 
£{a) = m, £{J3) — n a.s the lattice points of the region 

(2.21) %n,n = I {Xl, ...,Xm,yi,.-.,yn)e M>o+" ^ II ^' = II % h 

and for fixed g we may view the double Hurwitz numbers Hg(a,P) as defining a 
function 

/„ „„x / /IN I Aut all Aut /3I ^ , ^, 

(2.22) {a, /3) ^ J ^i-^ ^Hg{a, /3) 

on this set of lattice points. Goulden, Jackson and Vakil [22l Theorem 2.1] proved 
that there exists a hyperplane arrangement in ]R™+" and a collection of polynomials 
Pg, c{xi, . • . , Xm, yi, ■ ■ ■ , y-a) indexed by the chambers c of this arrangement such that 



the function (2.22) is given by 



(2.23) (a,/3)^Pg,,(ai,...,a„,/3i,...,/3„) 

for all (a, j3) e c. They used this piecewise polynomiality property to motivate and 
support a conjectural analogue of the ELSV formula for double Hurwitz numbers 
p2l Conjecture 3.5]. 

Following [22] , the piecewise polynomial structure of the double Hurwitz numbers 
was investigated by a number of authors [71 [Ml [571 [3S] . In particular, Johnson [33] 
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exhibits an elegant approach to piecewise polynomiality based on the formahsni of 
the infinite wedge space |51| . This representation-theoretic approach was recently 
extended by Shadrin, Spitz and Zvonkine [SF, where piecewise polynomiality was 
linked to the theory of shifted symmetric functions in a very structured way. We will 
use this viewpoint to elucidate the piecewise polynomial structure of the monotone 
double Hurwitz numbers. 

Let J-{Y) denote the algebra of functions Y — > C on Young's lattice, and consider 
the subalgebra A* of T freely generated by the functions 



(2.24) pUX)^J2 



i=l 



X.-^+l) -(-^ + 1 



k> 1. 



The algebra A* is known as the algebra of shifted symmetric functions, and the 
generators p^(A) are the shifted power-sum symmetric functions. Shifted symmetric 
functions /(A) are polynomial functions in the parts Ai, A2, . . . of the input partition 
which become symmetric after the change of variables A^ := Xi — i + ^. 

There are several ways in which the introduction of the algebra A* may be 
motivated. One reason for the ubiquity of this algebra in the representation theory 
of the symmetric groups is the Kerov-Olshanski Theorem ^J , which asserts that 
the central characters 

(2.25) /„(A) . Jit,; 

are a linear basis of A* . Another motivation is the deep analogy between random 
partitions and random Hermitian matrices [35], in which the shifted power-sums 
p^(A) play the same role as the moments tr(^'^) of a matrix. Finally, shifted 
symmetric functions provide a framework in which the precise relationship between 
the Hurwitz theory and the Gromov-Witten theory of P^ may be described [5T] . 

Consider the transform T : T{Y) — > J^(Y x Y) sending functions on partitions 
to functions on pairs of partitions defined by 

(2.26) T{f){a,P):^^^ ^ XaX^/(A). 

11 J llPi |A| = |a| 

This definition assumes that \a\ — |/3|; if this is not the case, i.e. if |a| > |/3| 
or vice versa, complete the smaller partition by adding an appropriate number 
of I's. For fixed m,n we may view pairs of partitions (a,/3) with \a\ = |/3| and 
£{a) = m,£(a) — n as lattice points of the region CH™,™, as above. For each pair of 
proper subsets / C [m] , J C [n] , consider the hyperplane 

(2.27) Wi,j = I {xi,...,Xm,yi,---,yn) G ^rn,n -^^^^^Vj 

This hyperplane arrangement is called the resonance arrangement in |34| I58j . A 
chamber of the resonance arrangement is a connected component c of £Hm^„\ U/cfml JcM ^i,J- 



The following general piecewise polynomiality property of the transform (2.26) re- 
stricted to the algebra of shifted symmetric functions is proved in [35] • 
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Theorem 2.2. To each triple (/, m, n) consisting of a shifted symmetric function 
f and a pair of positive integers m,n, there corresponds a collection of polynomials 
Pf ^ci^ii • ■ • 5 2;™, j/i, . . . , y„) indexed by the chambers of the resonance arrangement 
such that 

T{f){a,l3) = P/,c(ai,...,a,„,^i,...,/3„) 
for all {a, /3) G c. 

The proof of the piecewise polynomiality of double Hurwitz numbers follows 



directly from Theorem 2.2 Indeed, a straightforward argument as in [58] verifies 
that when (a,/3) € c, the transitivity condition in the definition of the double 
Hurwitz numbers is automatically satisfied. Thus double Hurwitz numbers are 
given by the character formula 

(2.28) L^^«("'^)-H^W E xix'.ihiW 

II- II- I'^l- |A| = |a| 

for all (a, /3) € c, where r = 2g — 2 + £{a) + £{(3) and /2(A) is the central character 
of the conjugacy class of transpositions in the irreducible representation of the 
symmetric group labelled by A. But this central character is a shifted symmetric 
function: 

(2.29) /2(A) -ip; (A). 
We thus have 



|Aut Q||Aut/3| ^ ^ /I 



(2.30) J ^, ^H,{a,^)^T{-p;] {a,P) 



for all (a,/3) G c, with r = 2g — 2 + £{a) + £{f3), so that piecewise polynomiality 



follows immediately from Theorem 2.2 



The above argument applies verbatim to the monotone double Hurwitz numbers, 
which are given by the character formula 

(2.31) ^7/^(a,/?) = ^i§l J2 xl^hMx) 

for all (a, jS) G c, with r = 2^ — 2 + £{a) + £(/?), once we verify that the function 
hr{X) :— hr{A\) is shifted symmetric. This fact in turn follows an alternative 
characterization, due to Kerov, of the algebra A* in terms of the contents rather 
than the parts of partitions (see [53], Proposition 2.4] for a proof). 

Theorem 2.3. The algebra A* coincides with the algebra of functions on partitions 
generated by p\{X) — |A| and the functions 

PkW--^PkiAx), fc>l, 

where Pki^-x) denotes the usual power-sum, symmetric function pj. G A specialized 
on the content alphabet A\ of A. 
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Indeed, one may check directly that ■^p'^W = PiW = hi{X), so that the char- 
acter formula for the usual double Hurwitz numbers is given by the T-transform 
of htir\ while the character formula of the monotone double Hurwitz numbers is 
given by the T-transform of hr- 
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